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analyzed. We find indications that the chiral symmetry breaking occurs for an inhomogeneous 
background as in the LOFF phase in condensed matter. 
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I. INTRODUCTION 



The idea of noncommutative space-time coordinates in physics dates back to the 1940's 
Recently, due to the discovery of Seiberg and Witten i2| of a map (SW map) that 
relates noncommutative to commutative gauge theories, there has been an increasing interest 
in studying the impact of noncommutativity on fundamental as well as phenomenological 
issues |3|J. 

Moreover, the idea of noncommutative space coordinates has been applied in condensed 
matter and in particular to the theory of electrons in a magnetic field projected to the lowest 
Landau level and to the quantized Hall effect |^. 

Another interesting feature of noncommutative field theories, also related with condensed 
matter, is that noncommutativity could represent a tool to describe the transition between 
ordered and disordered phases with inhomogeneous order parameters. 

In particular, the phase structure of theory has been recently discussed 
y, and, in 0, 1^ [l^ [ill, strong indications for a phase transition to a non- uniform 

stripe phase, due to noncommutativity, have been given. 

Originally the transition to an inhomogeneous phase has been considered in fermionic sys- 
tem to build a new non-uniform superconducting state in condensed matter (LOFF phase) 



12| . The interesting result is that the inhomogeneous phase can be more stable than the 
homogeneous BCS state with many relevant phenomenological consequences. This phe- 
nomenon has also been reconsidered in the analysis of the QCD phase structure and it has 
been proposed that, at large density, the QCD ground state is a color crystalline supercon- 
ductor [ij] that could be found in the core of a pulsar (for a recent review see jl^). 

In this paper we investigate if a noncommutative field theoretical model for interacting 
fermions shows a transition to a inhomogeneous phase where the order parameter, i.e., the 
fermionic condensate, is not constant in space-time. 

We consider the generalization of the cutoff Gross-Neveu (GN) model to 3+1 non- 
commutative coordinates and, by using the formalism of the effective potential for composite 
operators, introduced by Cornwall, Jackiw and Tomboulis (CJT), in the Hartree-Fock 
approximation we find that, due to noncommutativity i.e., 

there are indications for a transition to a non-uniform chiral symmetry breaking state similar 



to the LOFF one. 

The paper is organized as follows. In Section II we generalize the GN model to the 
noncommutative case and briefly review the CJT formalism ; Section III is devoted to 
a preliminary analysis of the occurrence of the transition to the non-uniform phase; the 
energy difference between the two phases is computed in Section IV and Section V contains 
some comments and the conclusions. 



II. NONCOMMUTATIVE GROSS NEVEU MODEL 

In this Section we shall summarize the formalism of the effective action for composite 
operators (see for details) and consider the simpler generalization of the GN model to 
the noncommutative case. 

For a fermionic field and for the composite operator, as < il){x)il){y) >, the CJT effective 
action T{G) is given by 

r(G) = -zTr hnSG-^ - iTiS-^G + T2{G) + ^Tr 1 , (2) 

where G{x, y) is the full connected propagator of the theory, 5* is the free massless propagator 

iS~\x-y) = -^5\x-y) (3) 

and T2{G) is given by all two particle irreducible vacuum graphs in the theory with propaga- 
tor set equal to G{x,y). The effective action is recovered by extremizing T{G) with respect 
to G and the Hartree-Fock approximation corresponds to retaining only the lowest order 
contribution in coupling constant to T2{G) (see 

We shall apply this formalism to evaluate the effective potential for the noncommutative 
generalization of the GN model which, in the commutative case, is defined by the chiral 
symmetric Lagrangian density: 

L{x) = iip^ip + g{ip-^)^ . (4) 

The canonical generalization of the model to the noncommutative case is obtained b 
substituting the standard product with the star (Moyal) product, defined as 



earned by 
1,.,4)Q 
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The effect of the star product on the Feynman rules of the theory is an additional momen- 
tum dependence in the interaction vertices for the "nonplanar" diagrams ( see while 
the "planar" diagrams have the same structure of the commutative theory. However, in 
the Hartree-Fock approximation of r2(G), the generalization in Eq. (0) does not introduce 
any "nonplanar" diagram due to the spin structure of the four fermion interactions and the 
corresponding calculation of the effective action is not different from the commutative GN 
case. 

Analogous to the noncommutative version of the 0{N) scalar model we can con- 
sider a more general expression for the noncommutative four fermion interactions which, in 
the planar limit, essentially reduces to the commutative GN model, but maintains genuine 
noncommutative contributions, i.e., nonplanar diagrams, also at lowest order in r2{G). 

The simplest generalization is obtained by considering the Lagrangian density 

In the standard case the addition of the second term is trivial since it reduces to a 
redefinition of the coupling g and to add a chemical potential contribution which disappears 
in the infinite volume limit. However, in the noncommutative case, it gives to r2(G), in the 
Hartree-Fock approximation, also a nonplanar term which introduces the noncommutative 
effects. In momentum space T2{G) turns out to be 




T2iG) = g[TiGip)TiGik){l + e'^'^P) - 2TrG'(p)G'(A;)] , (7) 

where k A p = k^O^^p^j and the traces are over all the quantum numbers. To obtain the 
previous expression for r2(G'), it has been assumed that the full fermion propagator G{x,y) 
is a translational invariant quantity. We shall comment on this point in the following section. 

The breaking of the chiral symmetry requires that the solution of the equation which 
minimizes the energy, 

is such that {G^p),^^} ^ 0. 
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It is impossible to study the transition to the new phase with the most general class 
of propagators G and we shall limit ourselves to a Rayleigh-Ritz variational approach 
where, however, a meaningful ansatz for G requires at least some physical indications on its 
asymptotic behaviors. 

First of all, let us remember that in the planar approximation, i.e., oo, where the 

noncommutative effects essentially disappear [Ij], the generalization proposed in Eq. © 
gives an analogous result to the standard GN model . 

n this case the translational invariant full propagator can be conveniently parametrized 
as l3 

where m is a constant which is determined by the minimum equation of the effective potential 

, [ d^p m 
m = Ag — -Vt^^ ^ • (10) 

However, for finite OA"^, it easy to check that the ansatz in Eq. is inconsistent with 
the minimum condition. Indeed, by inserting in Eq. Q the expression of T2{G) given in 
Eq. (jZj), and by using the parametrization in Eq. (jH)), the minimum equation for the mass 
turns out as (in Euclidean momenta) 



^ = 4^?/7^^^(l + 2e^^^^) (11) 



and the solution m = constant 7^ is ruled out by genuine noncommutative effects. 

Therefore, we first improve the previous ansatz in Eq. by introducing the following, 
translational invariant, parametrization of the full propagator 

^^"^'^^"V (27r)4p2_M(p2)2^ ' ^^^^ 

where the explicit dependence on the momentum has been introduced in the parametric 
function M(p^). 

Then, the minimum equation for M{p^) is (again in Euclidean momenta) 

2, r d^k M(k^) r d^k Mik"^) , , 



5 



To complete the Rayleigh-Ritz variational ansatz for the propagator and to evaluate the 
effective potential, one needs to know at least the asymptotic behaviors of the solution M(p^) 
for large and small (Euclidean) momenta. 

In Eq. (jl3p the p dependence is due to the second integral, since the first one is a constant 
for any function M(p^), which insures the convergence in the infrared region. 

However, the noncommutative term e^^^^ couples the infrared and ultraviolet asymptotic 
behaviors: due to the strong oscillating factor, for small p the integration region is dominated 
by large k and vice-versa. Then one has to proceed in a self-consistent way. One expects 
that, for large p, the noncommutative effects are negligible and a reasonable behavior is 

M(p2) ^ Me , (14) 



where Mq is a constant. Then, by Eq. (fT^ . one obtains 



m 



M(/)l,-«^8«/|^,^e-. (15) 
To simplify the calculations, the antisymmetric matrix 6^1, is assumed to be of the form 

9^" = 9!^^^ ^)®ld/2 (16) 
and 0, 0] the integration can be easily performed and it gives 



M(p^)Uo-|m,^, (17) 



which shows the leading behavior for small p, discussed in details in [6|, |9[ , due to the known 
IR/UV connection [17]. One can selfconsistently verify that, by inserting Eq. p7p in the 
gap equation Eq. ()15j) . the leading behavior of M(p^) for large p is a constant, as initially 
assumed. Then, a good ansatz for M(p^), which reproduces the asymptotic behaviors of the 
exact solution of the gap equation, turns out to be 

2g 1 



M{p') = Me 



1 + 



(18) 



Eq. ()12|) and Eq. p8|) represent the Rayleigh-Ritz variational parametrization of the full 
propagator and the constant parameter Mq has to be determined by minimizing the energy 
density. In this translational invariant case the relation between the energy density E and 
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the effective action Tti is well known llg and one has 



-E = — = I |21nll + "^^"y I -4 



J (27r)4 [ \^ p2 y p2 + M{p 



.2\2 



+ 



^7 (277)4 7 (27r)4[p2 + M(p2)2][p + M(F)2] V ^ ^ ' I 

where is the four-dimensional volume. 

As in the GN model, one finds the chiral symmetry breaking for gA"^ larger than some 
critical value {gA'^)^. The parameter Mq depends on the coupling constant and on OA"^, and, 
for — > oo, Mg — > m. 

However, the singular behavior of M{p'^) for small p suggests a possible non-uniform 
background and , as we shall discuss in the next section, the translational invariant propa- 
gator used so far should be considered as an approximation of a more deep dynamics. 



III. INDICATIONS FOR AN INHOMOGENEOUS CHIRAL SYMMETRY 
BREAKING PHASE 

As observed in the previous section, the leading behavior for small p of M{p'^) is ~ 
l/{9'^p'^). This signals ( despite of the translational invariant approximation) that the one 
particle irreducible (IPI) two point function is singular as p — and this physically amounts 
to a long range frustration: < ip{x)'ip{x) > oscillates in sign for large x 0,0,1^. Then the 
possible phase transition should be to an ordered inhomogeneous phase, where translational 
invariance is broken and the noncommutativity requires a nonuniform order parameter and 
a more general ansatz for the full propagator with respect to Eq. ()12|) . 

In the general case the order parameter is given by (a is the spinorial index) 

< i^{x)M^)^ >= / (^ / ^) (20) 

and it is a constant for the translational invariant case, i.e., G{p,q) = 6^{p — k)G{p). On 
the other hand, in the planar limit one has 

lim G{p,k) = 5\p~k)GM , (21) 

where Gpi{p) is the translational invariant solution of the planar theory in Eq. 0. 

Then, if one analyzes the problem for finite and large ^A^, where the noncommutative 
effects start ( let us remember that in cutoff unit p/A, k/A< 1 ), one can use the following 



approximation for G{p, k) 

Gip, k) ~ 5\p - k)GTi{p) + Fe{p, k) , (22) 

where Gti{p) is a translational invariant function which depends on 9h? and reduces to 
Gpiij)) for Oh? — > cx), and 

hm Fe{p,k) = . (23) 

Now, for X oo, the dominant contribution to < 'ip{x)aip{x)a > comes from the region 
p k i.e., 

< ^(a;)7A(x) >^ j ^[Gt/(p) + F,(p,p)] . (24) 

In other words, a translational invariant approximation mimics the right behavior for 

large x and, for large but finite 6', there is only a small deviation from the planar theory. 

Then, the results of the previous section give a good starting point to describe the fermionic 

condensate in these asymptotic regions, where one expects oscillating corrections to the 

constant background. This suggests the following form of the full non-translational invariant 

propagator to the order 1/^'^, in Euclidean momenta 

{-i> + M{p^))5\p-k) 
p^ + M(p2)2 + 

^ (5\p -k-P) + 5\p -k + P)) {P^A{p, k) + P^B{p, k)) , (25) 

where all quantities are expressed in cut-off units, the four-vector P = j and \P\ = 1, M{p^) 
is given by Eq. ()18p and A{p, k) and B{p, k) are, at this stage, generic functions. 

By replacing this propagator in Eq. in the Hartree-Fock approximation with inter- 
action given in Eq. ©, it turns out that, to order 0(1/6'''), 

V{G)=V{G)ti + P\/^T)nt . (26) 

where the non-translational invariant correction (Ar)7VT depends only on the function 
A{p, k), which we choose as 

Ah k) - + ^^^'^ (27) 



to preserve the spin structure of the translational invariant propagator 

After a straightforward calculation one obtains the following form of the correction 



V " J {2n)^J (27r)4[p2 + M(p2)2]p2 + M(/t2)2] 
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One should note that the remarkable factorization of the volume factor V, which follows 
from the ansatz in Eq. (j^^j) . despite of its non-translational invariance, implies that the 
right hand side of Eq. (j28|) is an energy density. 

For P ^ 0, due to the large oscillating factors in the integrands, the behavior of ( AF) ^rr 
is dominated by the integration regions of large p and k and it turns out that 



V 



(29) 



This result gives the indication of a transition to a nonuniform background related to 
the nontranslational invariant ansatz in Eq. (j2SI)- However, the CJT effective action has 
a clear physical interpretation as the energy density of the system, E, only for spacetime 
translational invariant propagators (see Eq. p9|) ). For static, but not space translational 
invariant systems, —T{G)\static = 't^Ti where r is the time interval and £t is the total energy 
[igI i. Therefore, the study of a possible phase transition to an inhomogeneous state, due to 
noncommutative effects, should, more correctly, be performed by taking 6*01 = 0, % 7^ 
with i,i = 1, 2, 3, and by using the time independent formalism. This is the subject of the 
next Section. 



IV. STATIC FORMALISM 



The static formalism has been developed in '16] only for the scalar fields and, in Appendix 
A, we extend it to the case of fermionic fields. The static propagator can be written as 



G{x,y}=G{0,x,y)=z 

J Ztt 



■ f + m 



(30) 



where F{x,y) = {p + m^)^/^(x, ^) and the two functions f{x,y) and m{x,y) describe the 
general time translational invariant solution of the gap equation (see Appendix A). 
The energy of the system £t is 

Tr(spin)(jr(x, X 

doj 



£t = -4m J d^x 
£x Tr(spin) 

static 1 



27r 



7o ujG{uj, X, x) — ■ VG(x, y) 



x=y 



(31) 



where F2 (G) [static corresponds to the last term in Eq. ((21), evaluated in the static limit. 



By following the same steps of the four dimensional calculation of Sections II and III, we 
initially consider, for the static propagator in the commutative case, the ansatz 

^(P) = ^0 (32) 

with constant m and = y/^~+mF, which gives the gap equation of the GN model in 
static limit 

If one considers the previous ansatz for G, in the noncommutative model in Eq. (jH)), the 
gap equation turns out to be 

which rules out a solution with constant m and requires a more general ansatz, where m 
has a parametric dependence on p, i.e., 

G{p) = ^'p'-'^-P + ^iP^ . (35) 

The gap equation is now 



with a selfconsistent asymptotic solution 

m{k) — > mo , k ^ CO , (37) 

where ttlq is a constant, the position 6'^^ = e^^^dk defines the vector 9 and x indicates the 
standard vector product. With the ansatz 

9 1 



m{k) = rriQ 



1 + 



(39) 



\k X 

and by following the same arguments given in Section III, it is straightforward to show that 
the energy of the system, evaluated by the nontranslational invariant ansatz for the static 
propagator 

Ge[p, k) = , — d%p -k) + 



2yp^ + ni{p)'^ 

]- (5\p- k-P)+ 5\p- k + P)) [P^A'ip, k) + P^B\p, k)) , (40) 
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where A' 



(p, k) is, analogous to the four dimensional case, 




(41) 



turns out to be lower than the noncommutative translational invariant case by terms of 



Therefore, the previous calculation gives a clear indication that the noncommutative 
effects are responsible for the occurrence of the chiral symmetry breaking in an inhomo- 
geneous phase, since the latter has always lower energy than the (translational invariant) 
homogeneous one. 

The qualitative agreement between the static calculation and the approach in Section 
III is expected due to the following points: i) we are considering only a slowly varying 
background with fluctuations amplitude suppressed by powers of 1/9; ii) the non-trivial 
factorization of the volume V in Eq. (PHjl makes possible a physically meaningful evaluation 
of the non-translational invariant correction to the energy density. 

In both cases (static and non-static), the energy density difference between these two 
phases is of order 0{g/6^), while the difference with respect to the planar theory is much 
larger. Then, for convenience, in Figure 1 we plot —r{G)/V, in the non-static calculation, 
for the planar theory and for the non-translational invariant case, and in Figure 2 it is shown 
that the noncommutative effects decrease the critical coupling constant with respect to the 
planar GN model. 

According to point i), one can easily evaluate the x dependence of the vacuum condensate 
to order 0(1/ 6^), which turns out to be 



where c is a constant and < ipip >ti is the constant order parameter evaluated in the 
translational invariant case. 

V. COMMENTS AND CONCLUSIONS 

Our computation of the CJT effective action shows that the chiral symmetry breaking 
occurs for an inhomogeneous phase, due to the noncommutative nature of the four fermion 
interactions considered in Eq. (jH)). The energy difference between the inhomogeneous and 



order g/O"^. 




(42) 
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the homogeneous phases, which both include the noncommutative corrections, is of order 
0{g/6'^). The order parameter has an oscillating x dependence of order Oi^g/O"^), superposed 
to the constant background of the translational invariant phase. 

These results are essentially based on the non-translational invariant ansatz for the full 
propagator in Eqs. (|25|) and (|27p. Let us notice that, in the commutative 1+1 dimensional 



GN model, non-translational invariant effects have been introduced in [2l| and a transition 
to the inhomogeneous crystal phase at non-zero chemical potential has been obtained. 

In y, ^ for the noncommutative scalar case, it has been observed that the boson con- 
densation does not occur in the mode = but there is a total depletion to k = Q where 
< (f){x) >oc cos{Qx). Analogously, in the fermionic case, our ansatz in Eq. (^3)) corresponds 
to Cooper pairs with a non-zero total momentum, as it happens in the LOFF phase in con- 
densed matter. The latter point provides an indication that the noncommutative cutoff field 
theory could be applied to describe the features of the transition to inhomogeneous phases. 



APPENDIX: A 

In this Appendix we derive the energy for time, but not space translational invariant 
fermionic systems with Lagrangian given in Eq. (jHI), in terms of the static propagator 
G{x, y) and we closely follow the procedure outlined in for the scalar theory. The total 
energy is related to the effective action, computed in the static limit 

£T{G{x,y))T = -T{G{x,y))Uatic • (A.l) 

The static limit of the effective action is obtained by taking the time translational invari- 
ant propagator G{x,y) = G{xo — yo,x,y) at equal time Xq = yo and by re-expressing F in 
terms of the static propagator, defined by the full propagator as G{x,y) = G{0,x,y). To 
obtain the form of G{x, y), we recall that the functional derivative of F with respect to the 
propagator is related, as explained in detail in to the bilocal source K{x, y) 

^^^'^'^^ - -U{xo-yo)K{x,y), (A.2) 



5G{xo - yo; x, y) 2 
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6I\ 
6G 



(A.3) 



and therefore, from the exphcit derivation of the effective action, one gets 

G~'^ix,y) = -i6{xo - yo)K{x,y) + S''^6^{x - y) + 

which shows that the most general form of G^^{x,y) is 

G'^{xo-yo,x,y) = S\xo - yo^o^ix - y) +iS{xo -yo)! ■ f{x,y) +im{x,y) , (A.4) 

where f{x, y) and m(x, y) are generic functions of the spatial coordinates and all the depen- 
dence of G~^{x, y) on the temporal coordinates is contained in the delta function 5{xq — y^) 
and its derivative . From Eq. ()A.4|) one gets the Fourier transform of G^^{xq — yo, x, y) with 
respect to the variable Xq — yo, 

/ + 00 . _, 
dxoe'-'^'"^G~^{xo, X, y) = -iuj-foS{x - y) + ■ f{x, y) + im(f, y) (A. 5) 
-oo . 

which can be functionally inverted 



G{uj, x,y) = i\ — (x, y) 



F2 



(A.6) 



where F{x, y) = {p + m^)^/^(x, y). 

Finally, the static propagator is obtained by integration 



G{x,y) = G{^,x,y) 



duj 
2^ 



F2 



[x,y) 



(7oF - 7 ■ / + m) 
2F 



[x,y) 



Incidentally, we note that the trace over the spin indices gives 



Tr 



(spin) 



G{x, y) 



/2m 



y) ■ 



(A.7) 



(A.8) 



We are now able to evaluate the effective action in the static limit and, for simplicity, we 
shall restrict the following calculation to a constant mass m{x,y) = m6{x — y). We start 
considering the first term in the general expression of F in Eq. namely — zTrLnG^^, 
where the trace refers both to spacetime and spin indices (we neglect in Eq. Q the logarithm 
of the free propagator S which gives a constant contribution to the effective action). With 
the help of Eq. ()A.5|) one gets [16 | 



d^x J — Tr(spin)Ln(G'~^(t^,f,f)) 



-2i 



oo 2tt 
°° duj 

o 27' 



dxQ J d^x 

2 J dxo J d^xF(x,x) = 4m J dxo J d^x[TT(^spm)G{ 



x,x) 



(A.9) 
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where we have used Eq. ()A.8|) to replace F in the last step. 

The second term to compute in Eq. Q is —iTiS^^G which yields, after integrating by 
parts, 



-i Tr(spin) J d'^x d'^y[p5'^{x - y)]G{xQ - yo; x, y) 



i Tr 



(spin) 



d X 



-^7o(-iu;)G(u;, f , x) - 7VG(0, f , y) 



4m / dxQ \ (i^x [Tr(spin)G'(a;, X 



iTr 



(spin) 



x=y 
'3 



dXQ / C? X 7 ■ VG'(x, 



(A.IO) 



a:=j/ 



Finally the term corresponding to r2 in the Hartree-Fock approximation are straightfor- 
wardly computed by replacing G{xq — yo, x, ^) 1^:0=^0 with G{x, y). By collecting the various 
contributions to the effective action, namely Eqs. ()A.9|) and ()A.10|) plus r2 in the Hartree- 
Fock approximation, we get the expression of the energy shown in Eq. (j3ip . 
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FIG. 1: The energy density E at [K^g/i-K^) = 2.5 for the planar theory (dashed hne) and for the 
noncommutative non-translational invariant case for 6A^ = 12 (a), 15 (b), 18 (c). 
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